Nonlocal response in thin-film waveguides: loss versus nonlocality and breaking of 

complementarity 
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We investigate the efTects of nonlocal response on the surface-plasmon polariton guiding proper- 
ties of the metal-insulator (MI), metal-insulator-metal (MIM), and insulator-metal-insulator (IMI) 
waveguides. The nonlocal effects are described by a linearized hydrodynamic model, which includes 
the Thomas-Fermi internal kinetic energy of the free electrons in the metal. We derive the non- 
local dispersion relations of the three waveguide structures taking into account also retardation 
and interband effects, and examine the delicate interplay between nonlocal response and absorption 
losses in the metal. We also show that nonlocality breaks the complementarity of the MIM and IMI 
waveguides found in the non-retarded limit. 



I. INTRODUCTION 



Guiding of light at mctal-diclcctric interfaces has at- 
tracted a lot of attention in recent years due to the sub- 
wavelength light confinement achievable by excitation of 
propagating surface-plasmon polariton (SPP) modesjii^ 
SPP guiding in a number of configurations is not limited 
by the diffraction limit, allowing for the manipulation 
and concentration of light on the nanoscalc)^ At the same 
time, stronger SPP mode confinement is typically asso- 
ciated with stronger mode absorption in metal, resulting 
in a trade-off between light confinement and propagation 
distancesi^ This trade-off can be tailored by consider- 
ing various waveguide structures, where especially waveg- 
uides based on thin metal films or narrow dielectric gaps 
between two metal surfaces have shown to provide a con- 
siderably better trade-offi^ Symmetric metal-insulator- 
metal (MIM) and insulator-metal-insulator (IMI) waveg- 
uides, see Fig. [Tl are the most fundamental of this class 
of waveguide structures, and provide a solid foundation 
for the understanding of more complex plasmonic waveg- 
uides. The key property of the IMI waveguide is its 
ability to support the so-called long-range SPP mode, 
which exhibits considerably low propagation loss. Fur- 
thermore, the MIM configuration forms the basis for the 
effective-index modeling (EIM) technique of more com- 
plex waveguiding structures, such as V-groove, slot and 
trench waveguides 

The MIM and IMI waveguides have been extensively 
studied experimentally^"— and theoretically<^ii^"— A key 
feature of any theoretical description of SPPs involves 
a suitable choice for the modeling of the response of 
free electrons of the metal. By far, the most common 
approach in the literature has been to apply the local- 
response approximation (LRA). The LRA solutions for 
the MIM and IMI structures were determined very early 
by Economoiiii and comprise two SPP modes being of 
even and odd symmetry with respect to the electric and 



magnetic fields, respectively. The properties of these 
modes are determined by their respective dispersion re- 
lations, i.e. by the relations between the frequency uj 
and the SPP propagation constant k, which are given by 
transcendental equations. In the nonretarded limit of the 
LRA, the surface modes of the MIM and IMI structures 
become identical)^ which is an interesting property that 
stems from Babinet's principle of complementary struc- 
tures 

However, issues with the approach of the LRA arise 
when either considering large values of k, where an un- 
physical limit is found for the frequency, or considering 
narrow insulator or metal widths {w < 10 nm), where 
singularities occur4i2£ In particular when investigating 
extremely narrow V-grooves, as recently realized experi- 
mentally;^ with techniques such as EIM, the inadequacy 
of the LRA manifests itself. Nonlocal response (or spa- 
tial dispersion) has been shown to remove this flawed 
behavior of the SPP modes in waveguiding structures 
such as single metal-dielectric interfaces^ infinite cylin- 
ders^ and more recently, hyperbolic metamaterials^^ 
conical tips^iS^ hybrid plasmonic waveguides^ wedges 
and V-groovesi^ 

Extensive theoretical work has been done on the fun- 
damental MIM and IMI waveguides, yet only a few stud- 
jgg22^-_3i jja,ve focused on nonlocal effects in these struc- 
tures. In this paper, we fill this gap by determining 
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FIG. 1. The three waveguide systems: metal-insulator (MI), 
metal-insulator-metal (MIM) and insulator-metal-insulator 
(IMI) along with the chosen coordinate system. 
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the dispersion relations of SPP modes of the IMI and 
MIM waveguides taking into account nonlocal response, 
retardation effects as well as interband transitions in the 
metals. We also revisit the simple metal-insulator (MI) 
waveguide structure. The nonlocal response is described 
by a linearized semiclassical hydrodynamic model,— 
which includes the quantum kinetics of the free-electron 
gas described by Thomas-Fermi theory. 

The derivations and corresponding results for the non- 
local retarded dispersion relations of the IMI, MIM and 
MI waveguides are presented in Sec. HH With the dis- 
persion relations for the three waveguides at hand, we 
examine in detail the interplay between losses and nonlo- 
cality in the metal by gradually increasing the absorption 
losses. The fundamental influence of losses versus nonlo- 
cality on the SPP dispersion has, to our knowledge, not 
yet been investigated. Furthermore, we compare modes 
of the IMI and MIM waveguides with and without retar- 
dation and nonlocality, and show that only in the non- 
retarded LRA do the modes of these two complementary 
waveguides become identical. Retardation and nonlocal- 
ity arc shown to break their complementarity. These top- 
ics arc discussed in Scc. lIIII Finally, Sec. |IV] concludes the 
paper. 



II. THEORY 
A. Nonlocal theory for thin-film systems 

To determine the modes of thin-film systems, we first 
outline the main equations for the electric and magnetic 
fields that must be solved. We then consider the class of 
guided solutions with transverse magnetic (TM) polar- 
ization. The boundary conditions for the metal-dielectric 
interfaces are also discussed. 

The free-electron gas of the metals comprising the thin- 
film waveguides is described by a nonlocal hydrodynamic 
equation of motioni^i^ An intuitive way of describing 
the effect of nonlocal response is that it serves to smear 
out the charges at the surface of the metal on the scale of 
the Thomas-Fermi screening length.— One of the key 
impacts of this charge smearing is the removal of field 
divergences that are known to occur in the LRAi^"— 
The hydrodynamic equation relating the current density 
J(r,aj) to the electric field E(r,a;) is given by^i^I 



spouse can be compactly written ai 



,22.34 



/3| 



V [V • J(r, u)] + J(r, w) = a(cj)E(r, u), (1) 



where (t(w) = i£oWp/(w + 17) is the Drude conductiv- 
ity, and = (3/5)i;p is the nonlocal parameter obtained 
from Thomas-Fermi theory, where Vp is the Fermi veloc- 
ity of the metal. By combining Eq. ([1]) with Maxwell's 
equations, the general equations describing the electric 
field E(r, uj) in a metal with hydrodynamic nonlocal re- 



f fc2 



V X E(r,w) = 0, 
)V-E(r,c.)-0, 



(2a) 
(2b) 



where fcm — fco\/em is the usual wave vector in the 



metal while fcni — \J + ^T'^ ~ Wp/eoo/Z^F is the ad- 
ditional longitudinal wave vector present in a nonlocal 
description of the metal. Here, fco = uj/c is the vac- 
uum wave vector, em = eoo(w) — uy^/{uP' + i^uj) is the 
local-response Drude permittivity including additional 
frequency-dependent polarization effects through eoo{^) 
not due to the free-electron plasma response. 

The electric field in the insulator regions with permit- 
tivity £d is described by the Helmholtz equation 



(v2 + fc2)E(r,c^) = 0, 



(3) 



where = k^yje^ is the wave vector in the insulator. 

Once the electric field has been determined, the mag- 
netic field H(r,a;) can be found from Faraday's law 



H(r,w) 



1 



ILDflQ 



-V X E(r,w), 



(4) 



and then the free-electron current density J(r,cj) in the 
metal can be found as 

J(r,w) = V X H(r,a;) + ia;eo£cx)(w)E(r, w). (5) 

Without loss of generality we set the propagation di- 
rection along the ^-axis and define the x-axis as perpen- 
dicular to the propagation plane, as in Fig. [TJ Then the 
electric and magnetic fields for TM polarization can be 
simplified to 



E(r,w) = [E,{x)e^+E,{x)e,]e'''\ 
B.(r, Lu) = Hy{x)eJ' 



(6a) 
(6b) 



where k is the SPP propagation constant. With the defi- 
nitions in Eq. ([6]), we can simplify the general expressions 
of Eqs. (I2][5]) to the following component form 



= 0, (7a) 

= 0, (7b) 

(7c) 
(7d) 



\dx'^ 




kE,{x) 


,dE,{x)' 
dx 


f 92 
\dx^ 




kE^x) 


J ,dE,{x)' 
dx 


Hy{x) =. 


1 


kE^{x) 


ox 


Jx{x) = 


-ikHy{x) + iuJSoeaoEx{x) 



which are to be solved in the metal regions, while in the 
insulator regions the governing equations are 



dx"^ 



E^x) = l^^^M. 
k dx 



(8b) 
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The magnetic fields in the insulator regions are also deter- 
mined using Eq. (|7c|l . For conveninence, we have defined 
a propagation constant normal to the interfaces in the 
respective regions given as 



Hj _ k fcj 



for j G {m, d, nl}. 



(9) 



With Eqs. ^ and solutions for the electric field, 
magnetic field and current density can be determined in 
the metal and insulator regions. At the metal-dielectric 
interfaces, we must connect the solutions using bound- 
ary conditions (BCs). Maxwell's BCs provide two of 
the three needed, namely the continuity of the tangential 
components of the electric and magnetic fields (Ez and 
Hy, respectively). In our treatment, we neglect effects 
due to electron spill-out and quantum tunneling, which 
unambiguously determines the third and additional BC 
to be the vanishing of the normal component of the free- 
electron current density ( J2.)i^iHi2id2 With this assump- 
tion we reduce our range of consideration to widths larger 
than 1 nm for the MIM waveguide42^ 



B. Single metal-insulator (MI) interface 

Before considering thin- film waveguides, it is instruc- 
tive to revisit the fundamental problem of SPPs propa- 
gating at a single MI interface. The MI problem with 
hydrodynamic nonlocal response in the metal has been 
solved by Boardman et at,-- in the simplest of cases 
where interband contributions and intraband damping 
were neglected. These results were recently generalized 
to include such contributions;^ however, without con- 
sidering the delicate interplay between the absorption 
losses and nonlocality in the metal, which we examine 
in Sec. |TTT1 

The retarded nonlocal dispersion relation for a single 
MI interface is exactly given as 



1 = Onl, 



(10) 



where S^i is an important nonlocal correction that will 
also appear below for the more complex thin-film waveg- 
uides, and is given as 



(11) 



We emphasize that when /^p ^ 0, the local- response dis- 
persion relation in Ref. 112 is retrieved since 6n\ —J- 0. We 
also note that the fc-solutions of the radical equation in 
Eq. (|T0|) can be obtained analytically by the standard 
method for solving radical equations through squaring, 
and can be represented in terms of the solutions of a 
third-order polynomial. 

For completeness, we also note that in the non-retarded 
limit c — > oo, the nonlocal correction 6^1 simplifies to 



(5"[ = lim 6ni 



k e„ 



Knl 



(12) 



and the dispersion relation for a single MI interface 
Eq. (fTO|) simplifies to 



1 



(13) 



The non-retarded local-response dispersion relation is re- 
trieved by letting in Eq. (fH)) . 



C. Metal-insulator-metal (MIM) 

The problem of determining the SPP modes of the 
MIM waveguide can, as in the LRA, be simplified by con- 
sidering the even and odd modes separately. The symme- 
try considerations apply to the electric field. The exact 
retarded nonlocal dispersion relation for the fundamen- 
tal, even mode is 



tanh 



(T) 



while for the odd mode we find 



(1 + Snl) 



coth 



£m ^d 



(1 + Snl) : 



(14a) 



(14b) 



where w is the width of insulator slab. These equations 
arc in agreement with the recent resultsi^ In the nonre- 
tardcd limit, Eq. (|14p simplifies to 



(15a) 
(15b) 



tanh('^)=-£i(l-H<5Sn, 
coth('^)=-£±(l + C)- 

^ / E in 



As previously mentioned, inclusion of nonlocal re- 
sponse regularizes the unphysical divergences encoun- 
tered in the LRA. This property is also preserved for 
the MIM waveguide, and we may see how it comes about 
by examining the limit w — >■ for the fundamental mode 
(which in the LRA produces a singularity). In this re- 
gard, we may neglect retardation effects and addition- 
ally simplify Eq. (|15a|) by using the small-x expansion 
tanh(x) ~ X. This yields 



2£d 



A. 



(16) 



where Amim is the nonlocal correction, which vanishes in 
the local-response limit f3-p — >■ 0. We emphasize that k 
stays finite even in the case of w = in contrast to the 
diverging local-response relatioi*^ given by Eq. (fTS]) with 

Amim = 0. 



D. Insulator-metal-insulator (IMI) 

The two SPP modes of the IMI waveguide can be clas- 
sified into even and odd modes, as in the case of the MIM 
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FIG. 2. Local and nonlocal complex dispersion relations of the SPP mode of the MI waveguide with ed = 1 and metal losses 
increasing from F = to F = 10~^. Local results are shown in red, while nonlocal are shown in green. Solid lines display the 
real part of the propagation constant, Re(A"), while the dashed lines display the imaginary part of the propagation constant, 
Im(_ft'). The value = 5 x 10"'' suitable for noble metals has been used. 



waveguide. However, while the symmetry characteriza- 
tion applied to the electric field for the MIM waveguide, 
here it is with respect to the magnetic fieldi^ The non- 
local modes for the IMI waveguide have previously been 
studied in the case of a lossless metal without interband 
contributions.— Here, we generalize these results to in- 
clude such contributions, which are important in realistic 
waveguides. The retarded nonlocal dispersion relation for 
the odd and even modes are 



coth 
tanh 



2 



— Onl coth 



£ d Km 



6n\ tanh 



2 

KnlW 



(17a) 
(17b) 



We note that, in contrary to the MIM waveguide, the 
odd mode Eq. (|17ap is in fact the fundamental mode. 
In the nonretarded limit, Eqs. p7a[) and (jl7bp simplify 

to 



' kw 
coth — 
2 

tanh I — 



- in rnr 



coth 



(18a) 
(18b) 



As in the case for the MIM waveguide, we can again 
examine the limit of w — >■ for the fundamental mode. 
Neglecting retardation effects and using the small-x ex- 
pansion coth(x) ~ 1/x, we find 

'A \ 21 

(19a) 

A,M,^T^J ^"^ (19b) 



2£d 


1 












2ed 






^nl^m V 





Here, Ami is the nonlocal correction for the IMI waveg- 
uide that vanishes for vanishing /3f. Unlike the MIM 
waveguide, the nonlocal correction does not regularize 
the diverging k when w = 0. Due to the confinement of 
the electron plasma in the IMI waveguide, as opposed to 
the MIM waveguide, the regularization of the dispersion 
likely requires inclusion of electron spill-out, which is not 
treated here, sec Sec. Ill Al Further elaboration on the 
comparison of Eq. (fT9|) with Eq. (fT6|) . in the context of 
complementarity, is done in Sec. IIIIBI 



III. RESULTS 

The dispersion relations introduced in Sec. |TT] are 
complex-valued transcendental equations of the implicit 
form F(uj, k) = 0, with the propagation constant in gen- 
eral being a complex number k = k' + ik" . Thus to de- 
termine the waveguide modes, solutions to the dispersion 
relations must be found in the complex fc-plane for each 
frequency, which in general is a non-trivial task. Fortu- 
nately, a robust and reliable numerical scheme suitable 
for determining the zeros in the complex plane, based on 
the Cauchy integral formula, has been previously devel- 
oped^i^ and is employed in this work. 

In the following, we focus on the free-electron prop- 
erties (i.e. £oo = 1) of the modes of the three different 
waveguides. This allows us to rescale the dispersion re- 
lations with normalized quantities, here introduced as 
Q. = w/ojp, K = kc/iOp, r = 7/wp, 77 = Py/c, and for 
the IMI and MIM waveguides, W = wuip/c. The nor- 
malized parameters F and 77 characterize the losses and 
the strength of nonlocality in the metals, respectively. 

This section is divided into two parts: Sec. IIII Al con- 
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FIG. 3. Nonlocal complex dispersion relations of the SPP mode of the MI waveguide with £d = 1 and metal losses increasing 
from r = 10~^ * to r = 10~^'^. Solid lines display the real part of the propagation constant, Ke{K), while the dashed 
lines display the imaginary part of the propagation constant, Im(_K'). The transition and mode-evolution from nonlocality to 
loss-dominated behavior is explored. The value rj — 5 x 10~^ suitable for noble metals has been used. 



corns the interplay between metal losses and nonloeal- 
ity. Here, we first study this interplay in the simple MI 
waveguide that does not contain any geometric length 
scales, whereafter we examine how nonlocal effects are 
enhanced in confined waveguides such as the MIM and 
IMI waveguides. Sec. IIII Bl deals with the breaking of 
complementarity in the MIM and IMI waveguides due to 
nonlocal response. 

A. Losses vs nonlocality 

MI waveguide: Due to the absence of length scales as- 
sociated with the geometry, the MI waveguide is an ideal 
system to study, when considering the interplay between 
losses and nonlocality of the SPP mode. An additional 
benefit of studying the MI structure is that it is not ob- 
scured by the effects of multiple interface reflections that 
is present in the MIM and IMI structures, such that only 
the intrinsic properties of free electrons affect the waveg- 
uiding properties. 

The interplay between losses and nonlocality in the MI 
waveguide is seen in Fig. [2l where we display the effect 
of increasing the metal losses on the local and nonlocal 
dispersion relations of the SPP mode, given by Eq. (fT0|) 
with 77 = 5 X 10~^. In the lossless case (F = 0), the lo- 
cal dispersion relation converges towards the well-known 
flsp = 1/V2ed hmit for large K values, while the non- 
local dispersion relation increases in frequency without 
bound, in agreement with earlier resultsi^ However, in 
the presence of very weak losses (F = 10"'^) the infinite 
K values at the frequency figp in the LRA are removed 
and the SPP mode bends back. This back-bending ef- 
fect is a well-known textbook result;^ which occurs for 



any positive value for F in the LRA. The extreme sen- 
sitivity to even minute losses in the LRA is due to the 
vanishing group velocity Vg = duj/dk at Osp41 In strik- 
ing contrast, the nonlocal SPP mode [i.e. Re(iir)] is ro- 
bust due to the finite group velocity Vg > /Sp. Conse- 
quently, no pronounced slow-light enhancement of weak 
losses takes place and the nonlocal SPP mode does not 
bend back until the losses of the system start to domi- 
nate. Although non-zero lm{K) is generated for the non- 
local SPP mode for F 7^ 0, the real part of the propaga- 
tion constant Ke{K) remains largely unaffected. It is also 
interesting to note that the behaviour of lm{K), which 
is related to the SPP propagation length ^spp through 
Zspp — l/[2Im(if)], changes drastically from F = 10"'^ to 
F = 10-2, For F = 10"^ the nonlocal SPP mode prop- 
agates longer than the local one in the frequency region 
n > risp, while the opposite result is seen for F = 10~^. 
At the same time Tie{K) for the nonlocal mode is un- 
changed and substantially larger than in the LRA, re- 
sulting in shorter wavelengths and thereby stronger con- 
finement of the SPP mode at the MI surface. Not until 
F = 10~^, which is significantly larger than the nonlocal 
parameter rj^ do the losses in the metal dominate over 
nonlocality and force the nonlocal SPP mode to bend 
back. At such losses, the local and nonlocal models re- 
sult in almost identical solutions. 

The transition of the nonlocal mode from being dom- 
inated primarily by nonlocality to being dominated by 
losses (i.e. F = 10"^ 10^^ in Fig. [2]) is investigated 
in more detail in Fig. [31 To explain the transition we 
must also consider the presence of the Brewster mode 
(for clarity not shown in Fig. and not only the SPP 
mode. In Fig. [3] we sec the merging of two separated 
modes, plotted as red and blue lines. For the lowest loss 
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of log(r) = —1.7, the red hne corresponds to the continu- 
ation of the Brewster mode to frequencies lower than D, — 
1 . which in the lossless case would be a forbidden region 
(i.e. only purely lossy solutions exist) 4^ The blue Hne 
represents the standard, low-loss, nonlocal SPP mode. 
As the losses increase [log(r) = —1.6 —1.575], the real 
parts of the dispersion of the Brewster mode and SPP 
mode begin to merge. At approximately log(r) — —1.55 
the mode-appearance has qualitatively changed, with the 
appearance of the usual well-known loss-dominated SPP 
mode (in blue), which is also present in LRA, as well as 
the emergence of a relatively flat band, nonlocal surface 
plasmon mode (in red) near the surface plasmon reso- 
nance fisp. We notice that the nonlocal flat band mode is 
significantly damped in comparison with the usual SPP 
mode, and that the damping increases drastically with 
increased material loss. In contrast, the usual SPP mode 
is not nearly so sensitive to the small change in material 
loss from log(r) = -1.55 to log(r) = -1.5. 

We now present a simple analysis to understand when 
the metal losses dominate nonlocal efi^ects in the MI 
waveguide. The back-bending occurs at the frequency 
r^sP) where the propagation constant is significantly 
larger than the free-space propagation constant. We can 
therefore justify to examine the simpler non- retarded dis- 
persion relation given by Eq. ()13|1 instead of the retarded 
dispersion relation [Eq. pH]) ]. From Eq. p^ . we see that 
nonlocality becomes negligible when |(5"f| ^ |1 + £m/ed|. 
Evaluating this condition at the SPP frequency with 
Ed = 1 (as in Fig. [2|) for small F leads to the simple 
condition for loss-dominated behavior 

F » 7?, (20) 

which is consistent with our numerical analysis. We point 
out that the loss parameter F is just one of several op- 
tions for introducing an imaginary part to the metal per- 
mittivity. An alternative approach to introducing losses 
is by simply adding a constant imaginary part to the 
lossless free-electron Drude model. In either case, the 
metal permittivity becomes complex-valued. To bridge 
these different approaches, we can relate the condition in 
Eq. (|20p to the imaginary part of the metal permittivity 
by noting that Im [em = i^sp)] 2V2T / (l + 2T'^) in 
which case Eq. PU]) can be rewritten as 

Im(£^) » ^. (21) 
c 

In noble metals the nonlocal parameter is of the order 
v-p/c « 10"'^ and the losses are of the order Im(eni) sa 10°, 
which in general means that metal losses largely domi- 
nate nonlocal effects in the SPP mode of the MI waveg- 
uide.— 

MIM and IMI waveguides: Figurc|4]displays the effect 
of increasing losses on nonlocality for the fundamental 
modes of the MIM (first row) and IMl waveguides (sec- 
ond row) given by Eqs. (|14ap and (|17a|) . respectively. The 
normalized width of the waveguides is set to W = 0.25 
which corresponds to a width of u; w 5 nm for Ag and 
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FIG. 4. Local and nonlocal complex dispersion relations of 
the fundamental mode of the MIM (first row) and IMI (sec- 
ond row) waveguides with W = 0.25, Ed = 1, J? = 5 x 10"'^ 
and metal losses increasing from r = to r = 10"\ Lo- 
cal results are shown in red, while nonlocal are shown in 
green. Solid lines display the real part of the propagation 
constant Re(Ar), while the dashed lines display the imaginary 
part of the propagation constant \m[K). The black lines rep- 
resent the approximate nonlocal dispersion relations given by 
Eqs. ([TeJ and ([TgJ for the MIM and IMI waveguides, respec- 
tively. 



Au (wp w 9 eV). Considering the MIM waveguide first, 
we see that in the lossless case nonlocal response within 
the hydrodynamic model predicts a blueshift compared 
to the LRA (for a fixed K). As the losses in the metal in- 
crease (F = 10~^), the local dispersion relation [Re(ii')] 
immediately bends back and the propagation length is 
significantly shorter than for the nonlocal case. Both 
of these effects are similar to those observed for the MI 
waveguide. When F = 10~^ the nonlocal dispersion rela- 
tion also bends back and the nonlocal propagation length 
becomes comparable to LRA, albeit for 57 < fisp nonlo- 
cal response gives rise to longer propagation lengths than 
in the LRA. Although the nonlocal dispersion relation 
bends back at these large losses, nonlocal response still 
reveals a blueshift and larger values of Re(i^) than in the 
LRA. 

The trend is very similar for the IMI waveguide (see 
second row of Fig. |4]). In fact, in the LRA the difference 
between the fundamental modes of the IMI and MIM 
waveguides is practically negligible. As for the nonlocal 
case, the biggest difference between the IMI and MIM 
waveguides is seen for F = 10^^, where nonlocal response 
shows a slight increase in the maximum values of both 
the Re(iir) and Im(ii') for the IMI waveguide. 

In Fig. m we have also examined the validity of the ap- 
proximate relations for the nonlocal fundamental modes 
of the MIM and IMI waveguides given by Eqs. (fTB|) and 
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([T9)) . respectively. They are plotted as black lines for the 
lossless case. We see that the approximate relations are 
in excellent agreement with the exact calculations when 
KW < 1. 

The important feature for both waveguides is that even 
for large losses (of order F = 10^^) the nonlocal and local 
dispersion relations arc different, in contrast to the MI 
waveguide. The nonlocal dispersion relations show larger 
values of Re(/^) than in the LRA for both waveguides. 
Thus, the limitations and undesired properties of metal 
losses are counteracted by nonlocality, which gives rise 
to a shorter wavelength of the SPP mode and thereby 
an increase of the mode confinement. These interesting 
features arise due to the multiple reflections present in 
the IMI/MIM waveguides, introducing a new length scale 
given by the scaled width of the slab W . The impor- 
tance of nonlocal effects increases with decreasing width 
(or in general, size);^ and it is clear from the nonlo- 
cal dispersion relations for the IMI and MIM waveguides 
that the strength of nonlocality is different in these two 
waveguide structures, as also observed in Fig. U) This 
difference arises due to the presence of confined nonlocal 
pressure waves, which are naturally only present in the 
IMI waveguide. 



B. Breaking of complementarity due to nonlocal 
response 

It is well known that the LRA dispersion relations for 
the SPP modes of the MIM and IMI waveguides arc iden- 
tical in the non- retarded limit j^ii^ which is also clear from 
comparing Eq. to Eq. (HI]) with 5™^ = 0. This prop- 
erty of identical surface modes in complementary waveg- 
uide structures, such as the MIM and IMI waveguides, is 
broken when retardation effects are included in the LRA, 
which become important for SPP propagation values K 
close to the light line = Here, we show explicitly 
that nonlocal response also breaks the symmetry by con- 
sidering the SPP modes of the MIM and IMI waveguides 
in the non- retarded limit, i.e. in the limit where K 3> Kt^. 
In the following, we divide the discussion of breaking of 
complementarity into two parts: one due to retardation 
effects alone in the LRA, and one solely due to nonlocal 
response in the non-retarded limit. For the latter, we 
consider the non-retarded limit to ensure that the break- 
ing of complementarity is due to nonlocal response rather 
than being attributed to retardation. 

Breaking of complementarity is illustrated in Fig. m 
In Fig. inja), we plot the SPP modes of the MIM and 
IMI waveguides only in the LRA, displaying the effect 
of retardation. First, we note, as already mentioned, 
that the SPP modes of the MIM and IMI modes in the 
non-retarded limit are completely identical and overlap 
in Fig. [HKa) (black lines). When retardation effects are 
included the MIM (green lines) and IMI (red lines) sur- 
face modes are no longer identical for K-M&hies close to 
the light line. The main consequence of properly taking 
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FIG. 5. Plots of both surface modes of the lossless IMI and 
MIM waveguides (a) with and without retardation in the LRA 
and (b) with and without nonlocal response (77 = 5 x 10~^) in 
the non-retarded limit. The light line is shown in grey. The 
widths of the waveguides are {&) W = 1 and (b) W = 0.25. 



retardation into account is that no guided modes exist 
above the light line (grey line). In Fig. [5^a), we clearly 
see that the retarded modes terminate at the light line, 
unlike the non-retarded modes. 

Figure mb) shows the nonlocal and local SPP modes 
of the MIM and IMI waveguides calculated in the non- 
retarded limit. We see clearly that nonlocal response dis- 
tinguishes between the MIM and IMI waveguide modes, 
for both of the two surface modes. This effect was ob- 
served earlierj^ but not elaborated on. As the propa- 
gation constant increases both of the nonlocal modes of 
both waveguides convergence towards the hydrodynamic 
nonlocal large- limit K = ^l/rj, as for the MI waveg- 
uide^^i^ and become indistinguishable. Finally, we also 
note the characteristic shift to higher frequencies of both 
nonlocal SPP modes compared to the LRA. 

The breaking of the complementarity property of the 
MIM and IMI waveguides due to nonlocal response can 
of course be understood from the fact that the dispersion 
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relations for the two waveguides are different even in the 
non-retarded hniit, as seen by comparison of Eq. (|15p 
with Eq. (fT5)l . We can quantify this difference for the 
fundamental mode of the two structures by considering 
the difference A/c = kuiu — kiuu where /cmim and kan are 
given by the approximate relations in Eqs. (|16|1 and 
respectively. Using a Pade approximation, we find to the 
lowest order in w that 

Afc~-^^. (22) 

From Eq. (22) we clearly see that in the absence of nonlo- 
cal response (Ajmi = 0), the difference between the funda- 
mental modes of the MIM and IMI waveguides vanishes. 
Additionally, we observe that the dispersion-difference 
depends strongly on the width. At very narrow widths, 
we therefore expect a strong breaking of complementar- 
ity. Thus, in the presence of nonlocality, a thin film of 
electron gas embedded in an insulator behaves qualita- 
tively different from a thin insulator gap embedded in an 
electron gas. 

In an intuitive, but simplified picture of nonlocality, 
one could be inclined to attribute the complementarity 
breaking to the nonlocal smearing of the induced surface 
charge. In the nonlocal hydrodynamic model the induced 
surface charge is smeared over a length scale compara- 
ble to the Thomas-Fermi screening length, leading to 
an effective width increase (decrease) for the MIM (IMI) 
waveguide. In this picture the dispersion relations of the 
nonlocal IMI and MIM SPP modes should then be be- 
low and above the local dispersion relations, respectively, 
which is not the case, see Fig. Hfb). In fact the nonlo- 
cal IMI and MIM SPP modes are always above the local 
results, discrediting the simple interpretation of nonlocal 
response as local response with effective size parameters. 

The complementarity breaking originates from the in- 
clusion of pressure waves in the description of a metal 
with nonlocal response. More precisely, the breaking is 
due to the confinement of these pressure waves in the IMI 
waveguide, which becomes more important for narrower 
widths. This confinement results in a significantly dif- 
ferent description of the IMI waveguide compared to the 
MIM waveguide, where the pressure waves are not con- 
fined. For this reason, nonlocal effects are also stronger 
in the IMI waveguide, as can be seen from the presence of 
only the nonlocal IMI correction in Eq. ([22|l. In contrast 
in the non-retarded LRA, the absence of both retarda- 
tion and the pressure waves leads to a faulty identical 



treatment of the MIM and IMI waveguides. 

IV. CONCLUSIONS 

The effects of nonlocal response, described by a lin- 
earized hydrodynamic model, on the waveguiding proper- 
ties of the MI, MIM, and IMI waveguides have been inves- 
tigated. The corresponding dispersion relations for the 
three waveguides have been derived, taking into account 
nonlocality, interband transitions, and retardation. The 
intriguing transition from nonlocal- to loss-dominated 
waveguiding behavior, which has not previously been 
studied extensively, was examined for the MI system, 
demonstrating that nonlocal response can counteract the 
effects of low metal losses. In the LRA the presence of 
even minute losses drastically alters the dispersion rela- 
tion of the SPP mode due the slow-light regime at the 
surface plasmon frequency i7sp- For larger losses, the ef- 
fects of nonlocality in the MI structure is less important, 
and the difference between local and nonlocal response 
becomes negligible. In general, for the MI structure, the 
impact of metal losses is much more pronounced than 
that of nonlocal effects, partially due to the high losses 
in metals and partially due to the absence of any geo- 
metric length scale in the MI structure. 

Conversely, for the MIM and IMI structures, the pres- 
ence of an additional length scale, given by the geometric 
width of the waveguide, yields a comparative boost to the 
effect of nonlocality vis-a-vis the effect of metal losses. In 
turn, the increased strength of nonlocality gives rise to 
larger propagation constants and thereby an increased 
plasmonic confinement of the SPP modes. Nonlocal ef- 
fects are shown to be slightly stronger in the IMI waveg- 
uide due to the presence of confined longitudinal pressure 
waves, which are absent in the MIM structure. 

Lastly, we also examined the complementarity prop- 
erty of the MIM and IMI waveguides in the context of 
Babinet's principle. In the non-retarded limit of the 
LRA, the waveguide modes of the MIM and IMI modes 
are known to be identical. When retardation is taken 
into account, this symmetry is broken. In addition, we 
have shown that in the nonretarded limit the symmetry 
is also broken by the inclusion of nonlocal effects due to 
the presence of nonlocal pressure waves. 
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